Effects of charge inhomogeneities on elementary excitations in La2 2:Sra;Cu04 
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Purely local experimental probes of many copper oxide superconductors show that their electronic 
states are inhomogeneous in real space. For example, scanning tuimeling spectroscopic (STS) imag- 
ing shows strong variations in real space, and according to nuclear quadrupole resonance (NQR) 
studies the charge distribution in the bulk varies on the nanoscale. However, the analysis of the 
experimental results utilizing spatially-averaged probes often ignores this fact. We have performed 
a detailed investigation of the doping-dependence of the energy and line width and position of the 
zone-boundary Cu-0 bond-stretching vibration in La2_a;Sri,Cu04 by inelastic neutron scattering. 
Both our new results as well as previously reported angle-dependent momentum widths of the elec- 
tronic spectral function detected by angle-resolved photoemission can be reproduced by including 
the same distribution of local environments extracted from the NQR analysis. 

PACS numbers: 63.20.-e, 74.25. Kc, 63.50.-x 



I. INTRODUCTION 

It is by now recognized that charge inhomogeneity is 
an important aspect of copper oxide superconductors. 
Undoped copper oxygen planes are Mott insulators due 
to a strong on-site Coulomb interaction. They become 
metallic and exhibit high-temperature superconductiv- 
ity (HTSC), when they are doped. "'^ Different types of 
charge inhomogeneity can emerge as a result of dop- 
ing. Nanoscale inhomogeneity due to strong correla- 
tions between electrons, such as a stripe or checkerboard 
order-i^, received a lot of attention in recent years. ^"'^ 
Other sources of inhomogeneity are randomly distributed 
heterovalent substituents or extra oxygen atoms, which 
introduce doped holes or electrons into the copper oxy- 
gen planes;^!^ In most cases (with some notable ex- 
ceptions, such as Ortho Il-ordered YBa2Cu306.5), the 
dopants act as charged impurities. They form an in- 
homogeneous Coulomb potential impacting the conduc- 
tion electrons in the copper oxygen planes, resulting 
in an inhomogeneous charge distribution. Local probes 
such as scanning tunneling microscopy (STM) have pro- 
vided evidence for this inhomogeneity in cuprates such as 
Bi2Sr2CaCu208+. (BSCCO),^-!^ Ca2-.Na,Cu02Cl2,13 
and La2-xSra;Cu04 (LSCO).^'* Bulk-sensitive nuclear 
magnetic resonance (NMR) and nuclear quadropole res- 
onance (NQR) experiments confirmed that an inhomo- 
geneous charge distribution is not just a property of the 
surface, but is a feature of the bulk in LSCO J^"— 

Here we show how some doping-dependent features of 
the observed spectra of phonons and electronic quasipar- 
ticles in LSCO, can be naturally explained in terms of 
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FIG. 1: Doping dependence of the Cu-0 bond-stretching 
mode energy, a) Dispersion in the [100]-direction at differ- 
ent doping levels, b) Zone boundary half-breathing {h = 0.5) 
mode energy as a function of doping (x). 



an inhomogeneous charge distribution. For phonons, we 
base our analysis on new results of neutron scattering ex- 
periments presented here. We use previously published 
angle resolved photoemission (ARPES) data for the elec- 
tronic response. 

It is a long-standing observation that the energy of the 
zone boundary Cu-0 bond-stretching (half-breathing) 
phonon, q = (0.5,0,0), depends strongly on doping. It 
changes from about 81 meV at x = to about 67 meV 
at a: = 0.3.—"— Figure 1(b), which combines previous 
neutron scattering and x-ray scattering data with our 
new results for the frequency of the half-breathing mode, 
summarizes the available data. Since the phonon fre- 
quency depends strongly on the doping, especially in the 
underdoped region, an inhomogeneous doping should in- 
crease its observed line width. We have measured the 
line width of this phonon in the underdoped region to 
investigate this effect. 
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II. EXPERIMENTAL DETAILS 

Inelastic neutron scattering experiments were per- 
formed on large, high-quality single crystals of LSCO 
with X = 0, 0.05, 0.07, and 0.15 on the IT triple axis 
spectrometer at the ORPHEE reactor at the Laboratoire 
Leon Brillouin at Saclay, France. The monochromator 
and analyzer were the (220) reflection of Copper and the 
(002) reflection of pyrolytic graphite (PG), respectively. 
The measurements were performed at reciprocal lattice 
vectors Q = {A + h,0, 0) in the tetragonal notation. In 
order to measure the intrinsic phonon line width, it is 
important to properly account for the experimental reso- 
lution, which depends not only on the instrument config- 
uration, but also on the phonon dispersion in the vicinity 
of the measured point in Q-w space. We calculated the ef- 
fective energy-resolution for the bond-stretching phonon 
branch at each doping level based on the spectrometer 
resolution combined with doping-dependent dispersions 
calculated from the shell model, whose parameters were 
optimized to agree with previously published phonon dis- 
persion data along high symmetry directions. The cal- 
culated experimental resolutions at Q=(4.5, 0, 0) for 
x=0.0, x=0.05, x=0.07 and x=0.15 are 4.5, 3.6, 3.2 
and 3.4 meV, respectively. Electron-phonon coupling 
or anharmonicity should result in a Lorentzian intrinsic 
phonon line shape, but we found that a Gaussian gives 
better fits. Thus, the resolution function was convolved 
with a Gaussian function representing the phonon peak, 
with the Gaussian width (and frequency and amplitude) 
determined by fitting to a constant-Q scan. The back- 
ground was determined from the scans at nearby wave 
vectors as was done in Ref. ^2 



III. RESULTS AND DISCUSSION 

Figure 2(a-d) shows inelastic neutron scattering (INS) 
spectra of the half-breathing mode of LSCO from x = 
to X — 0.15. The frequency clearly softens with increas- 
ing doping, as previously reported. Red lines through 
the raw data indicate the fitted Gaussian peak, convolved 
with the resolution function, on top of a linear back- 
ground. In the cases of a; = 0.0 and x — 0.05, we have 
included a second Gaussian to account for a small peak 
at 87 meV that we believe corresponds to an artifact, ^■^ 
based on its lack of dependence on doping, although we 
do not have a definitive explanation for it. We will ignore 
it in the following discussion of the line width; including 
it in the evaluation of the line width would not make a 
qualitative difference to the argument, as discussed be- 
low. 

The fitted phonon line widths are plotted as a function 
of doping in Fig. 2(f). The filled symbols show the width 
of the half-breathing mode at low temperature. As one 
can see, the width makes a large jump from x = to x = 
0.05, and then gradually decreases with further doping. 
For comparison, the open diamonds show the fitted width 
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FIG. 2: Inelastic neutron scattering data of half breathing 
zone boundary bond-stretching modes in LSCO. a) a; = 0, b) 
X = 0.05, c) a; = 0.07, d) a; = 0.15 at 10 K, and e) x = 0.05 at 
500 K. Red lines on raw data are composed of a linear back- 
ground and two Gaussian peaks for (a) and (b) and a linear 
background and a Gaussian peak for (c), (d) and (e). Large 
Gaussian peaks on the bottom represent longitudinal bond- 
stretching modes; small Gaussian peaks at 87 meV represent 
the a suspect artifact (see text). Full widths of half majci- 
mum (FWHM) of the main peaks of x=0, x=0.05, x=0.07 
and x=0.15 are 4.5 meV, 10.5 meV, 8 meV and 5meV, re- 
spectively. Two Gaussian peaks of x=0.0 data and a small 
Gaussian peak of a; = 0.05 data are resolution-limited (4.5 
meV). Dotted blue line on a: = 0.0 raw data is composed of a 
linear background and a Gaussian peak with 6.5 meV FWHM. 
f) Doping dependence of the zone boundary (h = 0.5) and 
zone center {h — 0.0) longitudinal bond-stretching mode line 
widths after deconvolving the resolution function. 



of the zone-center mode. The absence of any significant 
variation with doping suggests that the results for h = 0.5 
are intrinsic and are not a result of difference in sample 
quality. Furthermore, different samples with the same 
doping give identical results. 

The most obvious explanation of phonon broadening at 
the zone boundary is electron-phonon coupling. Accord- 
ing to the conventional theory of metals, the phonon line 
width due to electron-phonon coupling should be pro- 
portional to the density of states near Epr^^ which is 
roughly proportional to doping in LSCOi^"— The trend 
for the ft, = 0.5 mode in Fig. 2(f) clearly violates a sim- 
ple monotonic increase with doping. In particular, the 
X — 0.05 sample is insulating^ and should not have 
more conduction electrons near the Fermi surface than 
the X = 0.07 and 0.15 samples. Thus, we can rule out 
electron-phonon coupling as the dominant cause of this 
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FIG. 3: Doping distributions for x = 0.05 (a and g), x = 0.074 
(b) and x = 0.15 (c) based on the NQR data.^^ Inelastic 
neutron scattering data from the half-breathing mode for x = 
0.05 (d), a; = 0.07 (e) and a; = 0.15 (f) at 10 K and a; = 0.05 
at 500 K (h). Dotted red curves and solid blue curves over 
data points show results calculated from the inhomogeneous 
doping effect at high temperature and at 10 K, respectively. 
Red curve in (h) shows the simulated results at 500 K. 



effect. 

Another possible explanation of the broadening of the 
low-temperature line shape is the effect of the tilt of the 
CuOg octahedra responsible for the transition from the 
high-temperature tetragonal to the low-temperature or- 
thorhombic phase<^ We can rule out this scenario be- 
cause octahedral tilt is largest for x = and decreases 
with doping. 

Here we propose a different mechanism to explain our 
results. As already mentioned, NQR and NMR experi- 
ments found that the doping level in nominally homo- 
geneous samples of the LSCO family is strongly spa- 
tially inhomogeneous J^ii^ Since the phonon frequency is 
known to depend on doping, as shown in Fig. 1, we pro- 
pose that the observed broadening results from the distri- 
bution of local environments observed by NQR. To test 
this idea we modeled phonon line shapes based on the 
distribution of dopings measured by NQRiS. combined 
with the doping dependence of phonon frequencies. (As 



the intrinsic phonon line- width is smaller than the resolu- 
tion, we have ignored it in the modeling.) The resulting 
line shapes were convolved with the resolution function. 
This analysis provides a very rigorous test of the above 
hypothesis, since it contains no adjustable parameters 
aside from the overall peak magnitude. 

One can imagine more rigorous approaches to calcu- 
late the phonon spectrum in an inhomogeneous medium 
such as calculations based on the coherent potential ap- 
proximation (CPA) which is usually used for calculat- 
ing phonon spectra in alloys, But it is not clear if this 
method would apply to phonons in solids with charge in- 
homogeneity as in LSCO, since the CPA was originally 
developed for random, mass disordered alloysj^ How- 
ever, we note that calculations of disorder effects in iso- 
topically disordered Ge using CPA gave no better re- 
sults (actually even slightly inferior results) than calcu- 
lations using large supercells.^'^ Somewhat surprisingly, 
the supercell approach reproduced the experimentally 
observed phonon frequencies, linewidths and intensities 
very well, in spite of the fact that this approach is based 
on harmonic lattice dynamics. Therefore we think that 
our ansatz for calculating linewidths in LSCO - which 
is equivalent to a simplified version of the supercell ap- 
proach - is quite adequate. 

Figure 3 shows that the calculated line shapes re- 
produce the experimental low temperature phonon line 
shapes at three doping levels reasonably well. We also 
measured this phonon in the x — 0.05 sample at 500 K. 
Heating to such temperatures increases anharmonicity 
and thus broadens the phonon. For example, this broad- 
ening has been observed already at 330 K near the zone 
center in similar samplesj2S However, the 500-K data ap- 
pear to be slightly narrower than at 10 K, which is consis- 
tent with decreasing charge inhomogeneity with increas- 
ing temperature as reported by NQR.^- An anharmonic 
contribution must be present, but it appears to be com- 
parable to the contribution from charge inhomogeneity. 

It is not completely clear why charge inhomogeneity 
decreases upon heating the sample. One possibility is 
that some electrons are trapped by the disordered chem- 
ical potential due to randomly distributed Sr cationSfi^ 
Then, the real space distribution of the trapped electrons 
should broaden with increasing temperature resulting in 
a smoother disorder potential. Charge inhomogeneity 
resulting from electronic correlations should also decrease 
at elevated temperatures due to the breakup of electronic 
self-organization by thermal fluctuations. 

Based on the above evidence we conclude that the 
phonon linewidth can be well explained by charge in- 
homogeneities measured by NQR. Intrinsic phonon life- 
time, which is governed by anharmonicity and electron- 
phonon coupling, is too long to measurably contribute to 
the linewidth in addition to the effect of inhomogeneous 
doping. 

If this explanation is correct, charge inhomogeneity 
should show up in other experiments involving spatially- 
averaged probes. In particular, electronic quasiparti- 
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FIG. 4: Simulated Fermi surface results for (a) x — 0.07 and 
(b) X = 0.15. (c) Doping dependent Fermi surface evolution 
extracted from the ARPES data.— (d) MDC widths at Fermi 
energy from experiment -~ and simulation from anti-nodes to 
nodes. 



cles measured by ARPES should also be affected by in- 
homogeneous doping, since the Fermi surface in LSCO 
has a strong and nontrivial doping dependence: The 
Fermi wave vector k/ is nearly doping-independent at 
the nodes, but rapidly decreases with increased doping 
near the antinodes due to the proximity to a Van Hove 
singularity (Fig. Ac)^^ Similarly, measurements of the 
momentum width of the spectral function obtained from 
momentum distribution curves (MDCs) show that the 
width at the Fermi level is substantially larger in the 
antinodal region compared to the nodal onei^i^ We 
find that combining the chemical-potential-distribution 
model with the observed doping dependence of k/ pro- 
vides a reasonable explanation for the k dependence of 
the momentum widths. 

Figure 4(c) shows the doping-dependent Fermi sur- 
faces extracted from the ARPES datai2S. We simulate the 
ARPES data by introducing a Lorentzian broadening in 
energy of AE — 70 meV FWHM, which translates into 
AE/vp, broadening in k, at each k-point and assuming 
a distribution of the electronic band structures weighted 
by the inhomogeneous doping distribution, as we had 
done for the zone boundary phonon. 70 meV FWHM 
is a reasonable value based on the measured quasiparti- 
cle peak width near the Fermi energy in LSCO»^ We use 
tight binding bands with first, second and third neigh- 
bor hopping parameters fitted to the ARPES data as 
the electronic band structures for each doping .■iS. Figure 
4(a) and (b) show the simulated results for the Fermi 
surfaces of x = 0.07 and x = 0.15, respectively. Momen- 
tum widths of the quasiparticles near the Fermi surface 
broaden on going from the nodal to the anti-nodal di- 
rections for both dopings, as observed in experiments42 



Figure 4(d) demonstrates a good quantitative agreement 
between the MDC widths calculated from our model and 
the datai^ for x = 0.07. The fine width of ARPES data 
at the Fermi energy has been, so far, usually interpreted 
as the self-energy effect of the quasiparticle due to impu- 
rity scattering. There is, however, no reason that the line 
width near anti-nodal direction is much broader than the 
line width near nodal direction due to impurity scatter- 
ing which is believed to be isotropic. Our new analysis, 
yet, provides a natural explanation about broader line 
width near anti-nodal direction. 

An alternative model is to apply a broadening of fixed 
width in momentum rather than energy; however, it was 
found that uniform momentum broadening gives a much 
poorer agreement with the data. The physical signifi- 
cance of this result is unclear. 

Broader ARPES line widths near the anti-nodal di- 
rections have been interpreted in terms of stronger 
scattering,— but we quantitatively reproduced the ex- 
perimental data with a large enhancement of the MDC 
widths towards the anti-nodal direction by the inhomo- 
geneous doping effect. Our simulation, which only in- 
cludes doping inhomogeneity, cannot reproduce the en- 
tire ARPES spectrum, especially in the anti-nodal region, 
because of the pseudogap and the incoherent spectrum at 
higher binding energies.— However, we can conclude that 
the long-standing observation of broader peaks near anti- 
nodes in LSCO is due, in large part, to inhomogeneous 
doping. We find further evidence for the mechanism pro- 
posed here from another cuprate superconductor. Quan- 
tum oscillation measurements showed that an overdoped 
Tl2Ba2Cu06+d (T12201) has a much more homogenous 
chemical potential,— because the dopants in T12201 are 
farther from the copper oxide plane than the dopants in 
LSCO and BSCCO.^ Recent ARPES data on this mate- 
rial show that the spectral functions near the anti-nodes 
are even sharper than those at the nodes4^ These re- 
sults are consistent with the conclusion that the broader 
lineshape near the anti-node than near the node is not 
universal in the cuprates but is a material dependent ef- 
fect resulting from doping inhomogeneities. 

Crystal momentum is generally a good quantum num- 
ber for the electronic quasiparticles, whereas the length 
scale of charge inhomogeneities in real space is only a 
few nanometers .i^iiSrJ^ Thus, it may be surprising that 
the spectra of phonons and of the electronic quasipar- 
ticles are well reproduced by a linear sum of spectral 
peaks for each doping weighted by the inhomogeneous 
doping distribution function. For the electronic quasi- 
particles, 70 meV FWHM intrinsic line width, combined 
with the measured Fermi velocity, translates into about 
3 nm mean free path for the nodal direction, which is 
almost the same as the length scale of doping variation 
from NQR^— and STMii^ It is even shorter for the antin- 
odal quasiparticles, because they have a smaller Fermi 
velocity. We can, therefore, think that electronic quasi- 
particles reside in each patch, and simply add all spectral 
peaks from corresponding patches with different dopings. 
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For the phonon, we do not need additional line width 
broadening to fit the experimental data as shown in 
Fig. 3. This, however, does not mean that the phonon 
has a much longer correlation length than the electronic 
quasiparticles. If the wavelength of any waves in an in- 
homogeneous medium is smaller than the length scale 
of the inhomogeneity of the medium, the waves tend 
to be localized inside a domain of the inhomogeneous 
medium4i^— The wavelength of zone boundary phonons 
is two unit cells, which is about 4 times shorter than 
the length scale of the inhomogeneity measured by STM 
and NQR^ii^ Thus the distribution of the zone bound- 
ary phonon frequencies in a medium with inhomogeneous 
force constants directly reflects the distribution of the 
force constants. 

To summarize, we have found evidence that elemen- 
tary excitations in LSCO such as phonons and elec- 
tronic quasi-particles, are strongly influenced by inho- 
mogeneous doping previously reported by purely local 
probes, NQR, NMR, and STM. Phonon and electronic 
spectra are well reproduced by just averaging over spec- 
tra from patches having various doping concentrations, 
even though the typical length scale of a patch is a few 
nanometers. Good agreement between the model based 
on inhomogeneous doping with both ARPES and neu- 
tron results despite almost no adjustable parameters, 
provides compelling evidence that inhomogeneous dop- 
ing must be included into any analysis of experimental 
results in the cuprates. Recently, iron pnictide supercon- 
dutors have also turned out to have chemical potential 
inhomogeneities due to a random dopant distribution^^ 
It is, therefore, necessary to take account of these charge 



inhomogeneities when interpreting measurements on iron 
pnictide superconductors as well. 

Inhomogeneity most likely originates from a random 
potential created by introduced dopants (in LSCO it 
would be Sr), rather than from the self-organization of 
conduction electrons, for example into stripes. It is, yet, 
essential to gain further insight into interactions of dy- 
namic stripes with phonons, which is at the focus of a 
different project. If stripes interact with phonons via 
a mechanism discussed in Ref. f, then the effect of 
dynamic stripes should occur at a different wave vector, 
q — (0.25,0,0), not at the zone boundary. In fact we 
find that the huge phonon line widths reported near this 
wave vector in Ref. [s^ cannot be reproduced by inhomo- 
geneous doping alone for x > 0.07, and large additional 
broadening is required near optimal doping. These re- 
sults will be presented in a subsequent publication. 
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